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Abstract 

We investigate the possibility to describe nuclear matter in an approach constrained 
by the prominent features of quantum chromodynamics. We mapped the in-medium 
nucleon self-energies of a point coupling relativistic mean-field model on self-energies 
obtained in effective theories of QCD. More precisely, the contributions to the nu- 
cleon self-energy have been separated into the short range part, driven principally 
by the quark structure of the nucleon described in a quark-diquark picture, and the 
long range part, dictated by pion dynamics and determined using in-medium chiral 
perturbation theory. A saturation point, although unrealistic, is obtained without 
any free parameter. A realistic description of nuclear matter saturation properties 
has been obtained with the inclusion of a small phenomenological correction term 
to the short range part of the self-energy. 
PACS: 21.65. -Ff; ll.30.Rd; 12.39.Fe; 12.38.Bx 
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1 Introduction 



Effective field theories (EFT) based on hadronic degrees of freedom are at 
the basis of many successful models of nuclear structure. In particular, quan- 
tum hadrodynamics (QHD) models offer a field theoretical framework con- 
sistent with the symmetries of quantum chromodynamics (QCD). Effective 
Lagrangians, constructed with meson-nucleon or point-coupling interactions, 
have been applied successfully in the relativistic mean field approximation 
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(RMF) to describe nuclear matter and finite nuclei all over the periodic table 

m- 

The nuclear matter saturation arises in QHD models via a subtle cancella- 
tion between large Lorentz scalar and four-vector nucleon self-energies in the 
medium. There is evidence from nuclear spin-orbit splitting in finite nuclei 
that these self-energies should be of the order of hundreds of MeV. A recent 
many body calculation [3], using realistic nucleon-nucleon (NN) potentials, 
found that these large self-energies is a model independent fact, enforced by 
the structure of NN interaction. 

However, in QHD models, the Lagrangian parameters are determined in a 
purely phenomenological way to reproduce nuclear matter and finite nuclei 
properties. It would be desirable to get a clear connection of the nucleon self- 
energies to the underlying theory of strong interactions, QCD. On the one 
hand, this can be achieved, for example, with finite density QCD sum rules 
[4]. In this approach, the nucleon self-energy modifications in the medium are 
related to the modifications of the scalar quark condensate and quark density. 
However, these self-energies alone do not allow for a realistic description of 
saturation properties when taken into account in an effective hadronic model 
[5]. On the other hand, important progress have been made in describing the 
nucleon self-energies by using quark models, like for example the Nambu- Jona- 
Lasinio model [6j. This type of model is well suited for the study of short range 
effects in the mean-field approximation. 

Concerning the long-range part of the NN interaction, it is known long be- 
fore QCD that pion plays an important role. In the medium, the inclusion 
of pion exchange between nucleons in quark models, like the NJL one, leads 
to important technical difficulties and still seems out of scope. Another ap- 
proach is based on the concept of effective field theories [7j in which one has 
to write the most general Lagrangian compatible with symmetries of QCD. 
At low energy, pions and nucleons are the relevant degrees of freedom and 
the dynamics is controlled by the broken chiral symmetry. This is the basis 
of chiral perturbation theory [8] (ChPT), which has been successfully applied, 
for example, to construct NN potentials (see for examples [9|10fll| ). In this 
theory, the low-energy expansion of the Lagrangian is arranged in powers of 
derivatives and pion masses. The in-medium nucleon self-energies can in prin- 
ciple be determined with the in-medium ChPT as used first in [T2|. However, 
the determination of the contributions beyond the next-to-leading order in the 
medium is a really complicated task. Moreover, contact terms encoding the 
short distance dynamics arise at each order. 

In this work, we therefore investigate the possibility to obtain the nucleon 
self-energy in nuclear matter by using the complementary of these two ap- 
proaches. The one and two pion exchange contributions are determined at 
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next-to-leading order in the framework of in-medium ChPT. Instead of includ- 
ing the phenomenological contact terms, the short distance physics is described 
with a quark model of the nucleon. We use here a quark-diquark picture in a 
NJL model. 

This nucleon self-energy is then implemented in a nucleonic density dependent 
point coupling model in order to explore the possibility to describe saturation 
properties of nuclear matter. A saturation point is obtained, although at too 
low density, for some values of the quark mass without any free parameter. 
Reasonable saturation properties are obtained by adding a small correction to 
the quark-diquark picture. 

In section 2, we present the in-medium quark-diquark model of the nucleon 
in the NJL model. One and two pion exchange contributions to the nucleon 
self-energies are determined with in-medium ChPT in section 3. The nucle- 
onic point coupling model, the self-energies mapping procedure are detailed 
in section 4, together with the results and discussion. We conclude in section 
5. 



2 In-medium self-energies of the nucleon in a Nambu-Jona-Lasinio 
quark-diquark picture 

We derive in this section the in-medium modifications of the part of the nu- 
cleon self-energies associated to the internal quark structure of the nucleon. 
We use here a simple quark-diquark picture in NJL model. 

The Nambu-Jona-Lasinio [TSj model is a chirally invariant point coupling 
quark model, in which the dynamics is dictated by the dynamical breaking of 
chiral symmetry. The Lagrangian of two fiavor SU(2) NJL model reads in a 
general form: 

C^,i. = q[i^-mo\q + £}-'\ (1) 

where q is the flavour SU(2) quark field and mo the bare quark mass. £™* 
is the interaction Lagrangian, based on chirally invariant four Fermi interac- 
tions. Any four fermion interaction can be Fierz symmetrised and rewritten 
identically as a chirally symmetric linear combinations J2i giiq^iqY where gi 
are functions of the original couplings and Fj are matrices in Dirac, fiavor and 
color space. We consider the scalar, pseudoscalar and vector channels in or- 
der to reproduce the vacuum phenomenology of pion and omega mesons. The 
quark-quark interaction channels are similarly obtained by a decomposition 
of using Fierz transformation. For our purpose, we need only the scalar 
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diquark channel, which gives the leading contribution to the structure of the 
nucleon |T4]. We thus write the effective interaction Lagrangians in the qq and 
qq channels respectively: 



int,eff 



(qq) + {qil-,Tq) \- g^{q-i^qy 
q(3^'j''CT2f] Iq'C-^r^fi^l^q 



(2) 
(3) 



where C = 27270, and = y3/2A^ (A=3,5,7) are the color 3 matrices. 
The couplings and g^^, together with mg, will be chosen to reproduce both 
the vacuum pion and vector masses and the pion decay constant. The scalar 
diquark coupling Qsd is in principle related, via Fierz transformations, to the 
mesonic channel couplings. However, the form of this relationship depends on 
the choice of and is therefore not unique (see for example |l5]). As usual, 
we consequently have chosen the Qsd value to reproduce the vacuum nucleon 
mass, independently from the other couplings. 



2.1 Nucleon and mesons in vacuum 



The Dirac equation for a quark in mean-field approximation is given by: 



mo - 25(^70 (qi^q) + 2g^ (qq) q = 



(4) 



which defines a dynamical constituent-quark mass: 



m = mo - 2g^ {qq) 



(5) 



generated by a strong scalar interaction of the quark with the QCD vacuum. 
In the gap equation (Eq.(l5l)), the quark condensate {qq) can be written as: 



{qq) 



IrSik), 



where here Tr denotes traces over color, fiavor and spin. In Eq. 
resents the quark propagator defined as: 



(6) 

S{k) rep- 



S{k) 



m + ie 



(7) 
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The quark condensate is divergent due to the loop integral and requires an 
appropriate regularization procedure. As many authors [T6|T7] . we introduce a 
three-momentum cut-off A which will have the least impact on medium parts 
of the regularized integrals, in particular at zero temperatureflT]. In fact, since 
the model is non renormalizable, the cut-off A is just an additional parameter. 

As usual, the mesons are obtained by solving the Bethe-Salpeter (BS) equation 
in the quark-antiquark channels. Since this is the standard procedure, we only 
recall the principal steps. First, we define the quark-antiquark polarization 
operator in the M = vr, o", channel by: 



nM(g') = -^ [ T^Tr [TMS{k + q/2)TMS{k - q/2)] , (8) 

where the vertex = i'j^T^, il and stand for respectively the pion, 
sigma and omega mesons. Note that, in the vector channel, the Lorentz struc- 
ture of the polarization operator is n^'^(g^) = (^—g'^" + n^(g^). The me- 
son masses, ttim, and meson-quark-quark coupling constants, gMqq, are then 
determined by the pole structure of the T-matrix, i.e. by the conditions: 



l-KMnMiq^ = ml,) = 0, (9) 
'dUM{q^ " 



dq^ 

where Km = 2g.„, 2g.,, and 2g^ respectively for the vr, a and u mesons. 



(10) 



The nucleon is described as a quark-diquark bound state by solving the Fad- 
deev equation. The quark-quark bubble graph in scalar channel reads: 



n..(g) = 6z| ^Trz, [^,S{kh,S{q - k)], (11) 

where Tro denotes a trace in Dirac space. The BS equation in the scalar 
diquark channel is then given by: 



1 - 2(?,rfn,,(g) = 0. (12) 

The relativistic Faddeev equation is reduced to an effective BS equation for a 
diquark and a quark interacting via quark exchange when restricting the qq 
channels to the scalar one. We take the static approximation to the Faddeev 
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equation [14], where the exchanged quark propagator is taken as —i/m. The 
static approximation has been found to reproduce reasonably well the exact 
Faddeev result [18]. The Faddeev equation then reduces to a series of quark- 
diquark bubble graphs and the solution for the t-matrix in the color singlet 
channel is: 



with n7v(p) the quark-diquark "polarization": 



(27r) 



(14) 



in which r(g) is the diquark propagator 



4:igsd 



1 - 2gsdnsd{q) 



4:gsd + 



G 



sd 



(15) 



where the second equality results from the pole approximation of the propaga- 
tor used in the numerical calculation. We have verified that this is a very good 
approximation of the exact diquark propagator in the range of momentum 
considered here. In Eq. (fTSll . nisd is the diquark mass, and Gsd is defined as: 



(16) 



Note that for the determination of the polarizations (Eqs. (I8|), (fTTl) and (fT4l) ). 
we have used the same regularization procedure as for the quark condensate. 



2.2 Nudeon self- energy modification in the medium 



As outlined by Ripka in |19| . the in-medium modification of the quark con- 
densate is roughly the same in quark or nuclear matter in NJL model, as 
long as the vacuum constituent quark mass is not too high ( m ~ 400MeV). 
In first approximation, we thus consider the in-medium modification of the 
quark-diquark state in quark matter. The masses and couplings at finite den- 
sity are denoted with star superscript (as for example m* for the constituent 
quark mass at finite density). When working at finite quark density, one has 
to modify the quark propagator Eq. ((71): 
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• by making the replacements fc^j ^ = — 2g^ {<iln<i) and m m* 

• by adding to the quark propagator Eq. ((71) the medium insertion part 

-n^:^±^6{ko-E,)e{kF-\k\), (17) 

where kp is the quark Fermi momentum, and Ek = + ^g^p with p = 
(g7°g) the quark density and El = k^ = Vk^ + m*'^. 

The insertion of the medium part in each quark propagator generates finite 
integrals which can be performed almost analytically. The polarizations Ilj 
(j = M, sd, N) in the Bethe-Salpeter and Faddeev equations thus depend 
explicitly on the quark Fermi momentum kp- This explicit dependence in kp 
will not be indicated below. 

For a given Fermi momentum kp-, the product S{k)S{q — k) in Eq.lfTTil for the 
diquark is replaced in the medium by S{k*)S{qsd — k*), where = g° — ig^iP 
and qsd = q. If the equation 1 — 2gsdn.sdi(l^) = is satisfied, then the solution 
of 1 - 2gsdnsd{qsd) = with q^d = {\Jff + mf^ = g° - Ag^p, q) is, in the rest 
frame of the diquark, g° = m*^ + Ag^^p. 

Considering the nucleon "polarization" equation ffT4l) . S{k)T{p — k) is replaced 
in the medium by S{k*)T{pd — k*) with p° = — QguiP and pd = p. If the 
equation 1 + -^I\.n{p) = is satisfied at |^ = Mat + S^, where is the scalar 
self-energy, then the solution of 1 + ■:^Yiis[{pd) = is ^^^^ = + S° . We thus 

get p° = ^p2 + (Mtv + + Qg^jp- The scalar self-energy Sj^ is obtained 
by the solution of the quark-diquark pole condition, at rest (p = 0), in quark 
matter. The vector self-energy is Ey = Qg^iP- 

Repeating this procedure for each value of kp^ we thus obtain S^(pb) and 
2.3 Numerical results 

For a given cut-off A, or equivalently a vacuum constituent quark mass m, 
the free parameters mo, g-^ and g^^ are chosen to reproduce vacuum pion mass, 
pion decay constant and omega meson mass = 135MeV, /^r = 92.4MeV, 
= 782MeV. The parameter gsd is finally taken to reproduce the free nucleon 
mass Mjq = 939MeV. As usual, one has to consider several values for the 
vacuum constituent quark mass m. 

The results for the in-medium nucleon self-energies are reported on Fig. [T] for 
three values of the constituent quark mass between m = 370MeV and m = 
430MeV. For m < 370MeV, the nucleon is not bound and for m > 440MeV 
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the scalar self-energy of the nucleon is too much attractive. The parameters 
of the models and the values of the nucleon self-energies at a baryonic density 
Po ~ 0.17fm~'^ are given in Tab. [U 



200 
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-200 
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Fig. 1. Part of the scalar and vector nucleon self-energies from quark-diquark model, 
for three values of the vacuum constituent quark mass: m = 370MeV (solid line), 
m = 400MeV (dashed line) and m = 430MeV (dotted line). 

Some comments are in order, even if these results should not be interpreted 
alone. First, the scalar self-energy is attractive and the vector one is repulsive 
as expected. However, one interesting point is that there is a strong asymmetry 
between the scalar and vector part of the self-energy, with a ratio S^/Sy ~ 
-2. 



m(MeV) 


A(MeV) 


mo(MeV) 


5^A2 




gsd/g-K 


S^(MeV) 


E^(MeV) 


370 


605 


5.9 


2.3 


2.0 


0.735 


-260 


130 


400 


592 


6.0 


2.4 


2.4 


0.768 


-347 


160 


430 


583 


6.0 


2.4 


3.3 


0.789 


-490 


230 



Table 1 



NJL model parameters for three values of the constituent quark mass. The nucleon 
self-energies in the last two columns are given at the baryonic density po = 0.17fm~^. 
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Second, the vector self-energy is exactly linear in the baryonic density. As can 
be seen in Fig. [T] and verified numerically, the scalar self-energy is also linear 
in the density to a good approximation. 

Finally, the self-energies highly depend on the vacuum quark mass value. There 
is approximately a factor 1/2 from m = 370MeV to m = 430MeV for the 
part of the scalar self-energy coming from the quark structure of the nucleon. 
However, we will show later that the impact of the quark mass value on the 
nuclear matter saturation properties is not as crucial as it could be expected 
from these results. 



3 In-medium nucleon self-energies from pion exchange contribu- 
tions 

We derive in this section the in-medium modifications of the nucleon self- 
energies coming from the long-range pion exchange determined in the frame- 
work of chiral perturbation theory (ChPT). 

3. 1 Elements of chiral perturbation theory 

ChPT is an EFT of pions and nucleons, constrained by QCD chiral symmetry. 
Numerous reviews on ChPT applied to NN interaction (see for example [20]) 
or nuclear matter (see for example [T2f2T] ) are available. We will recall here 
only the principal features needed for our discussion. 

Since we need only pion exchange contributions, we consider only the vrA^ 
interaction Lagrangian. Since the interactions of the Goldstone bosons (pions) 
must vanish at zero momentum transfer and in the chiral limit (m^ 0), 
the low energy expansion of the ChPT Lagrangian is arranged in powers of 
derivatives and pion masses: 

^^ttN = + ^fl + • • • (18) 

The lowest order irN Lagrangian, with one derivative or pion mass insertion, 
can be written: 

C^l = i^iip - + ^7^75^^)^., (19) 
where tp is the nucleon Dirac spinor, and: 
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D^ = d^ + T^, (20) 

r. = ^ {^d,^ + ^d,e) = ^r.(7r X d,7z) . . . (21) 

Uf. = I {i^d.i - id,e) = -^T.d.Tz + ... (22) 

U = e = l + jT.n-^y + ... (23) 



with f/ = the SU (2) matrix collecting pion fields. 

The contributions in ChPT are usually analyzed in terms of the powers of small 
momentum over the large scale {Q/A^Y, where Q stands for the characteristic 
momentum of the problem under consideration (nucleon or pion momentum, 
pion mass, . . . ) and ~ 1 GeV is the chiral symmetry breaking scale. The 
chiral counting rule consists in determining at which power u each Feynman 
diagram contributes. 

For the self-energy diagrams we consider, the chiral counting rule gives rise to 
the following hierarchy of contributions: 

• = 3 is the lowest order (denoted order 0{Q^)), given by the one pion 
exchange diagram (one loop) with vertex insertion from 

• 1/ = 4 order (denoted order 0{Q'^)) receives contribution from the iterated 
one pion exchange diagram (two loops) with vertex insertion from C^^]^; 

• u = 5 order (denoted order 0{Q^)) gets contributions also from the iterated 
one pion exchange diagram, and irreducible two pion exchange diagrams, 
with vertex insertion from C^^li; 

• at z/ = 6 order, two pion exchange diagrams with one vertex insertion from 
the second order Lagrangian (with two derivatives or pion mass inser- 
tions) and the other vertices from C^^j^, contribute; 



Considering the technical difficulties of in-medium calculations, we will restrict 
ourselves to one pion (Fig. [2|) and iterated one pion exchange diagrams (Figs. 
[3] and SI). We thus will determine contributions at order 0{Q^). Irreducible 
two pion exchange diagrams contribute for only a few MeV [21], and we will 
thus neglect these contributions in this first approach. 

Finally, the Heavy Baryon formulation of ChPT [22] is often used for NN in- 
teraction. It consists in using a l/M^r expanded Lagrangian instead of jC^-tv- 
However, we will use the equivalent method of the Munich group |11| . One 
starts with the relativistic Lagrangian C^^lf and writes down relativistic am- 
plitudes. Relativistic vertices and propagators are then expanded in 1/Mj^. 
This method leads to the same results obtained when starting with a l/M^ 
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expanded Lagrangian, but is more efficient in dealing with the calculations. 



3.2 In-medium self-energy diagrams 



For the in-medium self-energy diagrams, the relevant small scale of the prob- 
lem is the Fermi momentum kp. Thus, the chiral hierarchy exposed above 
suggests that we determine contributions at order 0{kp). 

Denoting k = {kP = y^k^ -f- Mff, k) the four vector of the on-shell in-going 
nucleon, the self-energy for each of the diagrams considered can be written as: 



S(fc) = Es{k) + 7oSy (fc) - 7-kE,(fc) (24) 

As argued in [231124] . it is known from relativistic Hartree-Fock and Dirac- 
Bruckner calculations that the real part of the spatial component |k|St, is much 
smaller than S^, Sy. Moreover, in the case of pion exchange with pseudo- 
vector coupling considered here, |kE^|/E5 is of order |k|/M7v, with |k| < kp, 
for all the diagrams considered. We will thus neglect the spatial component 
S^, of the self-energy. 

We determine the in-medium contributions arising from the self-energy dia- 
grams using the procedure described in [T21I2T] . We recall here the principal 
features of this method. The diagrammatic calculation in nuclear matter in- 
volves the in-medium nucleon propagator: 



Sj^{k) = + M 



N 



^k^-Ml + ie 



2TT6{k^ - M%)d{ka)e{kF - \k\) 



(25) 



the pion propagator: 



and the tt — NN vertex from C)^j^. 



GM = ^— (26) 

'(1) 



fj-Jl.r. (27) 

As usually done [T21I25] . the pion exchange self-energy diagrams are consid- 
ered as perturbations to the vacuum. Consequently, the masses and couplings 
involved should be taken at their vacuum physical value. 
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The nucleon propagator splits additively into the vacuum SJf^ and medium 
S"™*"^ insertion parts. The calculation is organized according to the number 
of medium insertions. Self-energy diagrams with no medium insertion part 
contribute to the physical mass of the nucleon Mat. Therefore, we consider 
only diagrams with one or more medium insertions. 

3.2.1 One pion exchange diagram 

The one pion exchange Hartree diagram is trivially zero. The one pion ex- 
change Fock diagram, with one medium insertion on the internal nucleon line 
is shown on Fig. [2l The double line on diagrams stands for a medium inser- 
tion part of the nucleon propagator, the simple line for a vacuum part and the 
dashed line stands for the pion propagator. 



One has to write the complete amplitude of this diagram. After performing 
the integration on the time component of internal four vector, one is left with 
a three dimensional integral. The integrand is expanded in l/M^v, keeping the 
two first orders. The remaining integrals can then be achieved. The resulting 
self-energies are functions of both |k| and kp. As usual [26], we will take its 
value for a nucleon on top of the Fermi sea with |k| = kp. We thus will quote 
here only the expressions of the self-energy contributions for |k| = kp. The 
contribution to the nucleon self-energy of the one pion exchange Fock diagram 
(denoted S^., Sy) is given by: 



Fig. 2. One pion exchange Fock diagram. 




(28) 



with: 




(29) 
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and: 



where f = kp/rriT,. The convergence of the I/Mat expansion can be verified at 
this stage. The ml/Mjf contribution of S^. indeed represents 5% of the first 
order contribution, for a value of kp ^ 250fm~^ close to the saturation one. 

Finally, is given by: 

3.2.2 Iterated one pion exchange diagram 

The iterated one pion exchange Hartree diagrams with one or two medium 
insertions are represented on Fig.O The diagram with three medium insertions 
brings a purely imaginary part to the amplitude. The diagram A contributes 
to the physical pion mass. The diagram B is of order 0{1/M^) and is sub- 
dominant. Finally, there are four remaining one pion iterated exchange Hartree 
diagrams which contribute, namely C,D,E,F. 

We focus here on the C diagram with one medium insertion. First, the internal 
part of this diagram is the so called planar box graph which needs to be treated 
carefully. This is done in detail in [TT]. The amplitude involves four poles in 
each half complex plane. The procedure consists in performing the internal 
integral by Cauchy's theorem, and then expanding the remaining integrands at 
first order in 1/Mn for each pole. The dominant contribution arises from one of 
the nucleon poles. Using the tools of dimensional regularization procedure, it is 
found that this diagram is not divergent and does not need any regularization 
as outlined in [TT] . 

The introduction of a linear cut-off in the iterated one pion exchange contri- 
bution, as done in [T2|2T] . can be considered as an effective parameterization 
of contact terms. Alternatively, Lutz et al. choose to introduce explicitly a 
NN contact term in the Lagrangian in [25]. Either the cut-off or NN contact 
term is a parameterization of the short range physics. One of the goals of this 
work is to try to get a description of this short range physics more directly 
from the quark structure of the nucleon, and thus, to avoid double counting, 
we have not introduced such contact terms or cut-off in the pionic part of 
the self-energies. The short range physics will be explicitly taken into account 
in our approach with the quark-diquark description of the nucleon. It is thus 
clear that the interpretation of our results should be made only when both the 
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Fig. 3. Iterated one pion exchange Hartree diagrams with one or two medium inser- 
tions. 



quark-diquark and pion exchange contributions to the self-energy are taken 
into account. 



The contributions to the nucleon scalar self-energy of the one pion iterated 
exchange Hartree diagrams C (denoted S|) and (D+E+F) (denoted E|) can 
finally be written in the form, for |k| = kp'- 



-Gtt 



9Afn^ 



M, 



N 



fll 








■3 


4 2I 




- ' In 


1 + 4^;2 










24w 




.2 


+ 3" . 



arctan 



(2^)}, 



(32) 
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1 dw 



- 2 In (l + 4t; 



1 + Av'^w'^ 



2 2 
W 



W 



— w 



In 



+ dw 

Jo J —vw 



di 



W + l' 
w — 1. 



/I pv ^2 

rfu7 / —di 
-1 JO f 



+ 



2(1 + (t2) 



-ln(l 



cr 



In 



— V 



(33) 



with cr = ,^ti'+yf2 + ^2(y;2 _ j^') ^ defined earlier. The contribution to 
the vector self-energy of the nucleon is strictly equal to the scalar one. 

The iterated one pion exchange Fock diagrams with one or two medium in- 
sertions are represented on Fig. [H The diagram with three medium insertions 
brings a purely imaginary part to the amplitude. The diagrams A, B con- 
tribute to the physical masses and couplings. Finally, there are four remaining 
one pion iterated exchange Fock diagrams which contribute, namely C,D,E,F. 




Fig. 4. Iterated one pion exchange Fock diagrams with one or two medium insertions. 



The contributions to the scalar part of the nucleon self-energy from the iterated 
one pion exchange Fock diagrams C (denoted S|) and (D+E+F) (denoted S|) 
can be written as, for |k| = kp'- 
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M. 



Jo V 



yAnUj I Jo v^"" 1 + 2^2 

(1 + 8f + 8^^) arctan^ - (1 + 4^ ) arctan(20] } , 



(34) 



Jo 



1 + — In z ^ 



4^^ 1 + (^ _ ^)2 

+1 dw d^^ (AvV - ln(l + Av^w^)) 
(-Av^z" + ln(l + Av^z^)) 

+ dw d^^Ml + a') - a'] (in ^^"^ 



1 , 



R [i + {v + 0^] [1 + (^^ - 0^] - ^^^^ 



(35) 



with a — + as defined earlier, and: 



G{x)=v{l+v'^ + x^) 
1 

Ax 



__L[i + (, + ,)2][i + (,_,)2]i,l±|±||!^ 



R=^{l^y2_ p)2 + 4^2(1 _ ^2)^ 



(36) 

(37) 
(38) 



The vector part of the nucleon self-energy is again equal to the scalar part. 



3.3 Numerical results 



The physical input parameters are Mat = 939MeV, = 135MeV, /tt ~ 
92.4MeV and qa — 1-3. The self-energy coming from pion exchange is com- 
pletely parameter free. Even if the expressions obtained above for the self- 
energies contain all powers of the parameter fc^/m^, the chiral counting rule 
indicates that the diagrams taken into account contribute mainly at order 
0{kp), 0{kjp) and 0{kp). It is thus natural to take the scalar and vector 
self-energies as polynomial fits up to kp-. 
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1/3 2/3 

J ,jPb_,jPb_ 



N 



PB 

My 



(39) 



where the c* are density independent parameters in the i = S,V scalar or 
vector channels, and ps = -^kp is the baryonic density. These parameters 
are reported in Tab. [2l This numerical fit perfectly reproduces the exact self- 
energies. 





S 


V 


Co 


40.5 


40.5 


Cl 


469.9 


468.8 


C2 


252.9 


283.3 



Table 2 

Coefficients of the self-energies Eq. (f39ll for the pionic contribution in scalar (S) and 
vector channels (V). 

The self-energies are shown on Fig. [5] as functions of the baryonic density. 
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Fig. 5. Scalar and vector nucleon self-energies from ChPT at fifth order in kp as 
functions of baryonic density. 

Some comments are in order. First, the scalar and vector self-energies have 
very close values. This is not surprising since the iterated one pion diagrams 
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are taken at the non relativistic order C( I/Mat), which implies that = Sy. 
The small difference between the self-energies arises from the first relativistic 
(9(1/M^) correction taken in the one pion Fock diagram. 

Second, the part of the self-energies coming from pion exchanges is approxi- 
mately T,s,v{po = 0.17fm~^) ~ 130MeV, which gives a repulsive single particle 
potential. Let us recall that, at this point, these results include only the pionic 
part and that the short range physics is missing. 

Finally, higher order contributions in ChPT could be important. These con- 
tributions would arise first as order 0{kp) contributions. 



4 Nuclear matter 

We need a framework to describe nuclear matter in which the self-energies 
determined from nucleon structure in the quark-diquark model and pion ex- 
change contributions will dictate the dynamics. It should be covariant and 
flexible enough to impose the dynamics as obtained in the quark-diquark and 
pion exchange calculations. 



4.1 Density dependent point coupling model 

A point-coupling nucleonic Lagrangian with density dependent couplings of- 
fers a framework compatible with the requirements described above. Such a 
framework, with a mapping procedure of the self-energies, has been already 
used in [26] by Finelli et al. We are here interested only in the description 
of symmetric and infinite nuclear matter in the relativistic mean-field (RMF) 
approximation and the Lagrangian reads: 

-lGsipB)ii'i^)' - lGv{pB)ii'l,i^)^ (40) 

where is the nucleon field, is the baryonic density operator with psu^ = 
ip'j^ip and = (1 — v) is the quadri-velocity of the nuclear fluid. 

The coupling strengths Gs and Gy are supposed to be functional of the bary- 
onic density. These functions will be determined by mapping the nucleon self- 
energies on those obtained from the quark-diquark and pion exchange calcu- 
lations. This Lagrangian has to be used in mean-field approximation, with 
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fluctuations beyond the mean-field encoded in the density dependences. The 
equation of motion for the nucleon reads: 



[^^{id^ - V^) - {Mm + S)] = 0, (41) 
with the scalar and vector self-energies of the nucleon S and given by 



S = Gs{pB){i^^), 

dG 



u 



(dG. 
dpB 



dp 



(42) 
(43) 

(44) 



B 



In the RMF approximation for infinite symmetric nuclear matter, i/ji/j and 
ip'-j^ip are taken at their expectation values and the spatial components of 
currents vanish. The scalar and vector self-energies read: 



S = Gs{pb)Ps, (45) 

H = Gv{pb)pb + S„ (46) 

where ps = {'4'ip) and pb = (^tV) ^ire the scalar and vector densities. S,. is 
the rearrangement term in the mean field: 



„ 1 dG Q r, 1 OGv o 
2 dpB 2 dpB 

It has to be noted that the rearrangement term S,. is essential to ensure the 
density dependent models to be thermodynamically consistent. 

The self-consistent nucleon mass is related to the scalar self-energy: 



M-;, = Mn + S, (48) 

with Mat the free nucleon mass. The spin-orbit splitting in finite nuclei are 
strongly related to the value of the effective nucleon mass at saturation density 
Psat- To accurately reproduce spin-orbit splitting, it should be at least reduced 
to the value M]^(psat)/^Af ~ 0.64 |27j. We will consider that an acceptable 
range for the effective nucleon mass would be M^^^psai) / = 0.65 ± 0.05. 
The scalar density is given as a function of the Fermi momentum kp: 
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M 



N 



(49) 



where Bp = yikp)"^ + (^iv)^- The nucleon density is pb = -^kp. The en- 
ergy density and pressure of nuclear matter can be obtained from the energy 
momentum tensor 



P = EfP 



-GspI + -^GvPb, 

- £kin + -^GspI + -GvPb + ^rPB, 



(50) 
(51) 



where the kinetic term reads: 



£kin = -[SEppB + M*^p,]. (52) 

Empirically, the binding energy per nucleon, defined as £b = £ / Pb — M^, 
should present a minimum at the saturation density psat- At this density, the 
pressure vanishes -P(psat) = 0, which gives a non trivial constraint. We have 
considered the empirical range for saturation density psat = 0.153±0.015fm~'^, 
with binding energy £B{Ps&t) = —15.5 ± l.OMeV. 

We will also need the nuclear matter incompressibility modulus defined as: 



K = 9pBT^. (53) 

OpB 

The value of -ft'(psat) is related to isoscalar giant monopole resonances in spher- 
ical nuclei. However, the empirical estimates are strongly model dependent 
[28|29] and we will consider that a value of -ft'(psat) = 250 ± 50MeV remains 
acceptable. 



4-2 Determination of the density dependence of the couplings 



The preceding point coupling model is the framework in which we want to in- 
clude the quark-diquark and pionic exchange contributions to the nucleon self- 
energies. To achieve that, the procedure consists in mapping the self-energies 
in order to deduce the density dependent couplings Gs{pb) and Gv{pb)- This 
is achieved by writing: 
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Gs{pB)p3 



^5 + ^5' 



Gv{pB)pB + ^r = A + ^l 



(54) 
(55) 



where S^v are the self-energies deduced from quark-diquark picture in Sec- 

"Lr is defined 



'S,V 

tion 2 and are pionic contributions derived in Section 3 



Eq. (l47l) ). The self-energies coming from quark-diquark are taken in linear ap- 
proximation: 



^5 ~ 



G'sPs 



(56) 



yO 



(57) 



We have verified that, using an accurate fit in powers oikp for the scalar self- 
energy from quark-diquark picture, does not change the results. The pionic 
contributions have been written in powers of p^. We therefore have, consider- 
ing that ps ~ Pb up to saturation density: 



Gs{pb)=G% + 



1/3 2/3 

4 + + 



M. 



N 



Ml 



Gy{pB)pB + ^=\Gl + 



Ml 



1/3 2/3 



M. 



N 



Ml 



pB- 



(58) 
(59) 



The equation on Gv{pb) is a differential equation. The solution is: 



Gvipi 



G' + 



Ml 



6 



V 



^^s^Pb 



1/3 



-C 



6'^' M 



+ 



N 



2/3' 



{cy- 



1 snPb 



-2) 



Ml 



(60) 



4-3 Parameter free results 



We should at this point note that, as argued in [6], the NJL model generates a 
vacuum Mexican hat potential which contains an attractive tadpole diagram 
working against saturation. Consistently with the NJL model, we have verified 
that the nuclear matter cannot saturate if one takes into account only the 
quark-diquark part of the self-energy in the mapping procedure described 
above, neglecting the pionic contributions. 
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We now come to the results when both quark-diquark and pion exchange 
contributions are taken into account in the mapping procedure. In a first 
step, we have applied the mapping procedure, using diff'erent couples of values 
{Gg, Gy) corresponding to constituent quark mass values in the range m = 
370 — 440MeV. It has to be noted that, for a given constituent quark mass 
value, there is no free parameter to fine tune on the saturation properties. The 
binding energy does not present any minimum for the lowest quark masses, 
and presents an unrealistic saturation point for m = 440MeV, with density 
Pb ~ 0.03fm~^ and Sb ~ — O.SMeV. While it is far from the empirical region, 
this result is encouraging since, with no free parameter adjusted to reproduce 
the saturation point, it was highly non trivial to obtain a minimum. It is 
known that this minimum results from a subtle cancellation between the large 
scalar and vector self-energies. 

It is clear that our model suffers of some approximations, particularly the 
quark-diquark description in the NJL model. For example, the static approx- 
imation can lead in the vacuum to almost a 10% deviation from the exact so- 
lution of the Faddeev equation. An improvement of the quark-diquark picture 
would lead to a modification of the linear part of the self-energy. We explore 
the possibility in the next section to improve the nuclear matter description 
by taking into account a phenomenological correction to the quark-diquark 
picture of the nucleon. 



4-4 Phenomenological correction to the quark-diquark part of the self-energy 



In this section, we explore the possibility to improve nuclear matter saturation 
by adding a phenomenological correction term linear in the density to the self- 
energy. 

SJ:s = 6GIps, (61) 
5Ev = SG''yPB, (62) 



in the scalar (S) and vector (V^) channels. Since we get all the contributions 
to the self-energy associated to pion exchange at kp order, this correction 
should be attributed to the part of the self-energy associated to the quark- 
diquark dynamics. In order to avoid too many free parameters, we consider the 
same correction for the scalar and vector self-energies, namely 6Gg = 6Gy = 
6G^. The equations (I58|60lj are modified by adding this contribution to the 
couplings. 

For a given constituent quark mass value, 6G'^ is adjusted, via a least squares 
procedure to reproduce the saturation properties which we recall below: 
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• the saturation density psat = 0.153 ± 0.015fm~ 

• the binding energy Ssipsat) = —15.5 ± IMeV 

• the nucleon effective mass M'^{psat)/M]y = 0.65 ± 0.05 

• the nuclear matter incompressibility modulus K{psat) = 250 ± 50MeV. 

The results for different quark mass values are listed in Table [3l For each of 
these parameterizations, the average error on the four observables is less than 
one error bar, which is a remarkable improvement, since it has been obtained 
with only one free parameter. The binding energy versus baryonic density is 
represented on Fig. [6] for three quark mass values. 



m(MeV) (5G°(fm2) /),„j(fm-3) 8B{Me\) M*^/Mn i^(MeV) 



380 


-2.9 


0.134 


-16.3 


0.81 


249 


390 


-2.8 


0.134 


-16.2 


0.82 


245 


400 


-2.1 


0.128 


-16.4 


0.77 


260 


410 


-1.7 


0.122 


-16.3 


0.75 


267 


420 


-1.4 


0.12 


-16.3 


0.73 


274 



Table 3 

5G^ (see text) values and saturation properties of nuclear matter for different quark 
mass values. 

It is clear from Fig. [6] that the binding energy is almost independent of the 
constituent quark mass value for densities around and below the saturation 
point. Although there is only one free parameter, the saturation properties 
are reasonable for any of the quark mass values considered. However, the 
saturation density is somewhat too low and the effective nucleon mass too 
high. 

The best fit is given by the m = 410MeV parameterization with an average er- 
ror of 0.7 error bar. It will be denoted below PCMl parameterization. Figure[7] 
represents the evolution of the self-energies as function of kp for the PCMl pa- 
rameterization at the different steps of the calculation: quark-diquark+ pionic 
exchange contributions alone and the final result. It is clear that the most im- 
portant part of the self-energies is generated by the quark-diquark and pionic 
exchange contributions, whereas the phenomenological contribution is weaker. 
Quantitatively, this can be confirmed with the ratio SG'^/Gg, which is of the 
order 0.14. We can conclude that it is possible to obtain realistic saturation 
properties with relatively small correction to the quark-diquark picture of the 
nucleon. 

For comparison, in [3], the nucleon self-energies have been extracted at tree 
level from a realistic NN potential based on ChPT at fourth order. The authors 
found that most of the self-energies were generated by the contact terms. 
Similarly, in our approach, the quark-diquark picture is also responsible for 
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Fig. 6. Binding energy per nucleon versus baryonic density for three constituent 
quark mass values: m = 380 MeV (solid line), m = 400 MeV (dashed line), m = 420 
MeV (dotted line). 

the most important part of the self-energies. 

Till now, the part of the nucleon self-energies from ChPT has been deter- 
mined up-to and including order kp. At three loop level of ChPT, three body 
interactions appear and would contribute first at order kp. However, higher 
order efi'ects coming from the quark-diquark structure of the nucleon could 
also contribute at kp order, thus the inclusion of a phenomenological kp con- 
tribution raises ambiguity in its interpretation. The Appendix A is dedicated 
to the study of the effects of the inclusion of a phenomenological kp term on 
saturation properties and a comparison of the results with a Dirac-Brueckner- 
Hartree-Fock calculation is reported in Appendix B. 



5 Conclusion 



We have investigated the possibility to describe infinite and symmetric nu- 
clear matter in an approach constrained by quantum chromodynamics. We 
have mapped the nucleon self-energies of a point coupling relativistic mean- 
field model on self-energies obtained in effective theories of QCD. We have 
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Fig. 7. Evolution of the self-energies as functions of hp at the different steps: 
quark-diquark and pionic exchange contributions alone (dashed line) and self-en- 
ergy with the y?p correction added (dotted line) for PCMl parameterization. 

determined the contributions to in-medium nucleon self-energy by separating 
the short range part, driven principally by the quark structure of the nucleon, 
from the long range part, dictated by pion dynamics. 

We have taken the nucleon structure into account in a simple quark-diquark 
picture in a Nambu-Jona-Lasinio model, which is chirally invariant and re- 
produces the spontaneous chiral symmetry breaking. The quark-diquark pic- 
ture generates large attractive scalar and repulsive vector self-energies in the 
medium, with an average ratio at saturation density of E^/Ey ~ —2. 

The long range part has been determined by including the one pion and it- 
erated one pion exchange described in the framework of chiral perturbation 
theory at order 0{k\). The self-energies obtained are found to be approxi- 
mately 130MeV at saturation for scalar and vector channels. 

A saturation point with too low density is obtained without any free param- 
eter to fine tune for a given constituent quark mass value. The description 
of nuclear matter saturation properties (saturation density, binding energy, 
effective nucleon mass and incompressibility modulus) is improved by intro- 
ducing a correction term to the self-energy, linear in the density, which could 
be interpreted as a correction to the quark-diquark picture approximations. 
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This term is found to be relatively small, and we get a reasonable description 
of nuclear saturation considering that, for a given quark mass value, we have 
only one free parameter. 

To summarize, we have obtained a reasonable description of saturation prop- 
erties of nuclear matter with one free parameter. The possibility to obtain a 
rather good description of saturation properties by correcting quark-diquark 
picture approximations with one weak phenomenological correction term lin- 
ear in the density clearly indicates that, by refining the quark-diquark descrip- 
tion of the nucleon, the model could be noticeably improved. 



A Higher order chiral perturbation theory contribution 

We examine in this appendix the influence on the saturation properties of a 
kp contribution to the self-energies. At three loop level of ChPT, three body 
interactions appear and would contribute first at order kp. To our knowledge, 
there are not yet explicit calculations of these contributions available. We thus 
adopt a phenomenological point of view, introducing in the r.h.s. of equations 
(15811590 a new term to the self-energy defined as: 




in the scalar S and vector V channels. The implicit conjecture made in the 
ChPT is that these terms should be sub-dominant. We indeed expect that 
these new contributions, adjusted to obtain the best fit of saturation prop- 
erties, will be small. However, these new terms can include also higher order 
effects coming from the quark-diquark structure of the nucleon. 

In order to fine tune the saturation properties, the single particle potential 
U S + V is the pertinent observable to consider. The correction to U is 
approximately driven by the sum SGg + 6Gy. Thus, it does not matter how 
the kp contribution is split into scalar and vector self-energies. One can then 
take only one free parameter, and will obtain equivalent results by imposing 
different relations between SGg and 6Gy. We thus assume that the correction 
to the vector self-energy is zero. 



We have thus only one more free parameter 6G\;. The results are listed in 
Table [All for different quark mass values. 




SEl = 5G\pBp: 

5T.y=5Gyp\, 



(Al) 
(A.2) 



5Gl = 



(A.3) 
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-1.3 


-8.5 


0.156 


-15.4 


0.67 


288 


420 


-1.0 


-8.6 


0.146 


-15.0 


0.66 


292 


430 


-0.8 


-7.5 


0.138 


-15.2 


0.65 


306 


Table A.l 



, 5G\ values (see text) and saturation properties of nuclear matter for different 
quark mass values. 



The average error on saturation properties is less than 0.6 error bar in all 
cases. The best fit is obtained for m = 410MeV with an average error of 0.2 
error bar on saturation properties. This is a remarkable improvement by a 
factor four with regard to PCMl parameterization. The m = 410MeV set will 
be denoted PCM2 parameterization. 

The leading order correction term 6G^ has been reduced between PCMl and 
PCM2 from -1.7fm2 to -l.Sfm^. The ratio 5G7G^ is now 0.11, and the 5Gl 
correction is relatively small with ^G^Psat/G^s' ~ 0.12. The total phenomeno- 
logical correction is approximately of the same magnitude as with only a linear 
correction term. Let us mention that, in [3], the nucleon self-energies have been 
extracted at tree level from a realistic nucleon-nucleon potential. In particu- 
lar, the Idaho NN potential, based on ChPT up to fourth order, gives at third 
order a total pionic contribution of approximately the same magnitude as SG^. 

We can conclude that a kp contribution to the self-energy, even weak, could 
lead to an improvement of saturation properties. In order to identify the origin 
of this contribution, it would be valuable to obtain an evaluation of three 
loop level contributions in ChPT, so as to separate effects coming from pion 
exchange and from quark structure of the nucleon. 



B Comparison with DBHF calculation from realistic NN potential 

We present here the binding energy per nucleon and the nucleon self-energies 
on Fig. IB. II and IB.2I for the PCMl and PCM2 parameterizations. The results 
from a Dirac-Brueckner-Hartree-Fock calculation based on BonnA potential 
[30] are also shown. 

The empirical saturation point region is also represented on Fig. IB.ll as a 
dashed area. It is clear that PCMl parameterization leads to a too low sat- 
uration density, while the PCM2 saturation point is approximately in the 
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Fig. B.l. Binding energy as functions of the baryonic density for PCMl and PCM2 
parameterizations and DBHF calculation based on BonnA potential (stars). The 
dashed area is the empirical region of the saturation point. 

centre of the empirical region. The DBHF calculation [30] gives a too high 
saturation density psat = 0.185fm^^, but it is, to our knowledge, the most 
elaborate DBHF calculation based on a realistic NN potential. As can be seen 
from Fig. IB. 11 the incompressibility modulus is lower for DBHF calculation 
(K ~ 240MeV) than for our models. 

As the three calculations give rather different saturation densities, we found 
more instructive to compare the self-energies plotted against the Fermi mo- 
mentum, normalised to the saturation Fermi momentum of each model. This 
is done on Fig. IB.2[ As can be seen, the PCMl parameterization gives weaker 
values for the self-energies, but the PCM2 parameterization leads to a good 
agreement with the DBHF BonnA calculation for densities in the range Pb = 

0.1 - i.epsat- 
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